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ABSTRACT 

Simple analytical expressions may be derived for sound 
velocities in cubic crystals by using lattice harmonics or 
functions which are invariant under the crystal symmetry operations. 
Thesd expressions are in good agreement with the exact results 
for typical crystals such as metallic iron and potassium 
fluoride . 
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The propagation of sound waves in single crystals is governed 
by a set of three simultaneous linear equations known as the 
Christoff el equations. [l-3 ] The characteristic equation 
relates the sound velocity v to the direction n of sound 
wave and to the elastic constants of the crystal. When the 
direction n lies along a crystal symmetry axis, simple solutions 
may be obtained for the sound- velocities . In many cases, these 
simple solutions may be readily reversed so that the measured 
sound velocities may be used to calculate the elastic constants. 
However, there are many situations where the solution of the 
characteristic equation is required for arbitrary crystallographic 
directions? for example, an argon crystal grown at low temperatures 
may be inaccessible except for sound velocity measurements. 

For these arbitrary directions, it is advantageous to have 
simple approximate solutions, preferrably analytical in form. 

Many experimental quantities, such as nuclear magnetic 
resonance (NMR) second moments and sound velocities, are invariant 
under the crystal symmetry operations, thus they may be expanded 
in terms of lattice harmonics of the identity representation. 

These lattice harmonics are specific linear combinations of 
spherical harmonics which are invariant under the crystal symmetry 
operations. [6-8J The linear expansions are very convenient 
for certain experimental quantities such as NMR second moments 
where it can be rigorously shown that the expansion is limited 
to only the zeroth, second and fourth order spherical harmonics. 

[ 9 ] However, such limitations may be absent for the sound 
velocities, thus the linear expansions many converge rather 
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slowly. In this paper, the simple analytical solutions for sound 
velocities along the major crystal symmetry axes have been 
extended to arbitrary directions for cubic crystals. However, 
we have not restricted ourselves to linear expansions alone. 
Instead, we have used simple functional forms which are chosen 
to fit the experimental data exactly along various directions 
with high crystal symmetry. These approximate solutions have 
been compared with the exact results, and the agreements are 
satisfactory. 

By choosing the reference or coordinate axes to be coincident 
with the edges of a cubic cull, then the cubic crystal is 
characterized by three elastic sti fness constants, c^ , c^ 2 
and c^» The Kelvin-Christoffel stiffnesses are given by [1-3] 

F ii = n i. 2 ( c n “ + c 44 (1) 

Fij = n i n j ^ c i 2 + C W for i ^ J < 2 ) 

where n^ , n 2 and n^ are the components of the unit vector n. 

The sound velocities v are obtained by solving the characteristic 
equation or determinant J"T^j - pv 2 . j = 0 where p is the 
density of the crystal. This equation may be simplified by 
introducing the dimensionless variables A =(p v -c^)/(c^-c^) 

and P = ( c i2 +c 44^ c ll -c 44^ so " tha " t eigenvalue A depends 
only on a single parameter 0 characteristic of the cubic 
crystal! 
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On expanding the determinant, we get the cubic equation 

A 3 -X 2 + p(i -<s 2 )a- q(i -p) 2 (i + 2(3) = 0. (4) 

22 22 22 222 
where P = iij n 2 + n g n^ + n^ n^ and Q = rij n g n^ are fourth 

and sixth order functions which are invariant under the crystal 

symmetry operations. These functions are equivalent to the lattice 

harmonics [ 6 , 7 ] after orthogonalization and normalization; for 

example, the fourth order lattice harmonic is given by 

-(525/16tt)*[p - (1/5)]. From (3) or (4) , the eigenvalues have 

simple forms when n is parallel to one of the three major crystal 

axes, the cubic edge [100] , the face diagonal [lioj, and the body 

diagonal Tin]. [1-3] 

n || [lOO] : P=0 , Q=0 , A l =1, A t =0,0. 

n (| [lio] i P=l/4, Q=0, A l =(1+/9)/2, A t = 0, (l-/3)/2. 

n || [ill] i P=l/3. Q=l/27, A L =(l+2p)/3, A T =(l-/*)/3. (1-/3)/^ . 

The largest eigenvalue A^ refers to the longitudinal wave wher as 

the other two eigenvalues refer to the two transverse or shear 

waves. It may also be seen from (3) or (4) that A T =0 for one of 

the transverse waves when n is perpendicular to one of the cube 

edges (Q=0) . Furthermore, simple results may also be obtained 

when n is perpendicular to one of the face diagonals. By defining 

\ = u(l-(3) , then (4) may be written as 

[l-£J 2 f(u 3 ~ u 2 +Pu-Q) + ^9 (-u 3 +Pu-2Q)] =0. (5) 

If u is independent of (3 , then both parentheses in (5) would be 

'<2 1 

zero. Solving u -u +Pu-Q=0 and -u^+Pu-2Q=0 for P and Q, we get 

P = 2u - 3u 2 , Q = u 2 (l-2u). (6) 

When n is perpendicular to one of the face diagonals, then n ^ =n £ 

2 2 

and n-j = 1 - 2n^ . Thus we have 
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P = 2^ 2 - yn^ 9 Q = n^U - 2n x 2 ) . (7) 

2 

Comparison with (6) indicates that u=n^ is one of the solutions. 

Thus one of the transverse wave eigenvalues is A T = (1-/9), 

For n^=n 2 , the characteristic equation is obtained by the substitution 
of (7) into (4) » 

A 3 - A 2 + Pd-/* 2 )^ 2 - 3n 1 4 )A- n 1 ^(l-2n J 2 ) (1-^) 2 (1+2/3) = 0. 

2 

Since one of the solutions is n^ ( 1 — /3) , this cubic equation may 

p 

be factorized into the product of A-n^d-^3) and the quadratic 
equation 

^ 2 - [l+n 1 2 (p-l)]A - (|3-l)(2/3 +l)n 1 2 (l - Zn^) = 0. (8) 

The other two eigenvalues may be obtained by solving (8). 

We will now extend these simple solutions to arbitrary 
directions of n. For the longitudinal wave, the simple function 
A L = 1 - 2(l-j3 )P (9) 

agrees with the eigenvalues for n along any one of the three 
major axes (100, 110, 111). The next order correction would be 
proportional to an invariant function up to the eighth order. 

This function would be a linear combination of P, Q and P . By 
imposing the condition that this function vanishes when n is 

r<j 

along one of the three major crystal -axes , the result is 
4P 2 -P-3Q. Thus we have the result 

\ L = 1 - 2(1-/3)P + K(4P 2 -P-3Q). (10) 

The constant K will be chosen so that A L agrees with the exact 
value when n is along the [112] crystal axis (this crystal axis 
is perpendicular to both a face diagonal and a body diagonal) , 
where we have n^ 2 = x\ 2 - 1/6, r\ 2 2/3, P = 1/4, Q = 1/54, and 
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4P 2 -P-3Q = -1/18. From (8), we get A L = [5 + ^+(9-6/3+33^ 2 )^]/l2. 
Substitution into (10) gives K = (3/2) [ 1 +5j3- (9-6^+33p 2 )^] . 

Thus we get 

X L = 1 - 2(1-/3)P + (3/2)[l+5^-(9-6j9+33p 2 )*J(4P 2 -P-3Q) (11) 

for the longitudinal waves. 

From (3) or (4) , the sum of the three eigenvalues is always 

unity, hence it is necessary to evaluate A T for only one of the 

transverse waves. We note that Arp=0 when Q=0 (n perpendicular 

to a cube edge) and that Arp = n^ (1-^3) when (n perpendicular 

to a face diagonal). Because of these restrictions, it may not 

be advantageous to expand Aj as linear combinations of lattice 

harmonics or invariant functions. For n perpendicular to a face 

2 2 2i 2 

diagonal, we get P -3Q = (n^ - 3n^ ) from (7), hence we have 

n^ 2 = P - (P 2 -3Q)^ and it follows that 

A t = (1-(9)[P-(P 2 -3Q)^ ] (12) 

This relationship is also exact for n perpendicular to a cube 

"V 

p 

edge ( \ T =0 when Q=0). We note that from (7), n^ may also take 

2 2 

other forms. For example, we may also write n^ = (2P -6Q)/(P-9Q). 
However, the expression A,p = (1-^) (2P 2 -6Q)/(P-9Q) will not 
vanish for Q=0. 

The anisotropy (or the variation with the direction of the 
sound wave) of the eigenvalues depends only on the dimensionless 
parameter ^ . For fi= 1, the three eigenvalues are A =1,0,0 for 
any direction* that is, there is no anisotropy for the sound 
velocities. Usually (although not always), we have f3< 1 for 
ionic crystals and |3>1 for metals; the extreme values are p =0.19 

- 5 - 


for SnTe and ^=6 for metallic sodium. [2, 10 ] Typical values 

are ^=2.1866 for metallic iron (c^^=23»30» 0^2=13*92 an ^ c 44 = i 1# 62 
11 2 

in units of 10 dynes/cm ) and =0.523 for potassium fluoride 
(c^1=6.58, c^ 2 = 1*^9 and c^=1.28). The exact eigenvalues for these 
crystals are shown in Fig. 1 and Fig. 2 as the solid lines for n 
in the (100), (110) and (111) crystal planes by solving (3) or (4) 
numerically. The points in the figures are the approximate 
eigenvalues obtained from the analytical expressions (11) and (12) 
for the longitudinal and transverse waves. For all orientations 
of the sound wave direction n, the agreement is quite good between 
the exact and the approximate eigenvalues. 

In conclusion, it may be seen that simple analytical expressions 
may be used with good accuracy for sound velocities in cubic 
crystals. These expressions are based on the lattice harmonics 
or functions which are invariant under the symmetry operations of 
the crystal. Whereas linear expansions can be used for the 
longitudinal waves, other functional forms may be necessary for 
the transverse waves because of the "cross-overs" of the two 
branches. These simple results may be useful for fitting 
experimental data of materials of unknown orientations crystallized 
at low temperatures or for the inversions of sound velocity data 
to obtain the elastic constants. [4,5] 
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FIGURE CAPTIONS 

Fig. 1. Eigenvalues A^ and A^, of the longitudinal and transverse 

waves for metallic iron ( ^=2.1866) when n is in the (100), 
(110) and (111) crystal planes. Solid lines are the exact 
values from (4) . The points are the approximate values 
from the analytical expressions (11) and (12). 

Fig. 2. Eigenvalues A^ and of the longitudinal and transverse 

waves for potassium fluoride ( ^= 0 . 523 ) when n is in the (100) , 
(110) and (111) crystal planes. Solid lines are the exact- 
values from (4) . The points are the approximate values 
from the analytical expressions (11) and (12). 
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